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Abstract 

Iterating Newton's method symbolically for the general quadratic ax 2 +bx + c yields a rational function 
A , the numerator and denominator of which are polynomials with highly composite coefficients. In 

particular, the coefficients have no prime factors greater than 2 n after n iterations even though most of 

the coefficients are much larger than 2 n . 

1 Introduction 

If / : C — » C is a diffcrcntiable function, then Newton's method, applied recursively to an initial value of zoj 
yields the sequence of values z\, Z2, ■ ■ ■ defined as 

z n+1 = N(z n ) = z n (1) 

f (Zn) 

which, in many cases, converges to a root of /. We restrict our attention to the general quadratic f(x) 
c'.r + bx + c for the purposes of this paper. Instead of using equation Q as a numerical method, we are 
interested in iterating symbolically by letting the initial guess be zq = x where x E C Doing so, we obtain 

_ P n {x) 

Qn{x) 

where P n (x) and Q n (x) are relatively prime polynomials in C[a;]. An interesting observation, pointed out 
to us by Jim Propp, is that the coefficients of these polynomials, though very large, have no prime factors 
> 2™. On his suggestion that this divisibility property might imply some combinatorial interpretation not 
immediately apparent in the formulation of equation Q , we have derived an explicit symbolic formula for 
each iterate. 

The aim of this paper is to show that the polynomials P n and Q n are given by the following explicit 
formulas: 

P n ( x ) = a 2 -V" - £ Yl (— 1) J ( . ) ( ■ )a k+j b 2 "- k - 2j - 2 c i+1 x k (2a) 



2"_1 2"-fe-l 



J 



Qn(x)=Y, e (-^uir V 1 )" L ji ' 2 L 1 ( 2b ) 

fe=0 3=0 V / V / 

The binomial coefficients immediately explain the maximum size of the prime divisors of the coefficients of 
these polynomials. Furthermore, equation J2J provides a good starting point for investigating any combina- 
torial meaning of P n and Q n , which is discussed at the end of this paper. 
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2 Fractional Linear Transformations 



Let C* denote the Riemann sphere C U {oo}. A fractional linear transformation is a map /3 : C* — > C* of 
the form 

ct + a 

where a,b, c, d € C and ad — be ^ 0. This function is a conformal map which is analytic everywhere except 
the pole at t = — = . For our purposes, we define a particular fractional linear transformation <^{t) as 

<p(r) = (3) 
r - r 2 

where n and r 2 are the two distinct roots of the quadratic polynomial f(x) — ax 2 + bx + c. This particular 
approach exploits the fact that Newton's method for quadratics (with distinct roots) is conjugate to z — > z 2 , 
with respect to tp, the fractional linear transformation which sends the roots of the quadratic to and 
oo. Rick Kenyon [2] was the first to point this out to us; expositions can be found in both Cayley PP and 
McMullen [3]. For distinct roots r% and r 2 , the fractional linear transformation <^ _1 (r) exists and the formula 
from equation Jy| can be expressed as 

ajy \ az n + bz n -\- c _-y 2 

N M = Z »- 2az n + b = <P {lf{Zn)) 

That is, the following diagram of maps is commutative: 

C - > C 



C* > C* 

(■) 2 

Thus, the formula for the nth iterate of Newton's method is 

^-l(W x )2 n ) = ^i(x-r 2 ) 2 " -r 2 (x-ri) 2 " = P n (x) 
(x - r 2 ) 2 " - {x - ri) 2 " Qn(») 

Using this formula, we shall prove the equations in J2J). 

3 Proof of Explicit Formula 

Theorem 1. Given a quadratic polynomial f(x) = ax 2 + bx + c, with b 2 — Aac ^ 0, define polynomials P n 
and Q n as follows. 

a 2 " -1 / „\ 

P»(x) = (ri(x-r 2 ) 2 -r 2 (x-ri) 2 ) (5a) 

r\ — r 2 \ J 

Qn(x) = ( (X - T 2 f - (X - nf) (5b) 

n - r 2 V / 
T/ien i/ie polynomials given in equation |5|) are egwa/ to i/iose in equation 0) which are reproduced below. 

2—2 2 —k— 2 /on\ /on _ i. _ _ o\ 

P n (x) ^a 2 "-^ 2 " - MM fc )( ■ )a fe +^ 2 ™- fe - 2 ^V +1 x fc 

fc=o j=o V / V ^ / 



2 ' ''" 1 / - 1 " > '2™ - k - j - 1 



fc=0 j=0 
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Before we begin the proof of the theorem, we will need the following lemma. 
Lemma 1. For all x, y S K, and all n £ N, n> 1, the following identity holds 



x n -y n 



2/)£(-l) 



i=0 



n — i — 1 



{x + y) n -^-\xy) 



(7) 



Proof of Lemma^ We proceed by induction. For n = 1,2 the equality is easily verified, so assume the 
identity is valid for all k < n. Let T(n) be the right-hand side of equation J7J. Since 



x n+1 - y n+1 = (x + y){x n - y n ) - xy(x n - x - y^ 1 ) 

we need only prove that 

T(n + 1) = (x + y)T(n) - xyT(n - 1) 
Combining the sums on the right side of equation JSJ by shifting indices, we have 

n-l 



(8) 



- y) 



n — 1 




n — i — 1 

i 



n — i — 1 

i - 1 



(x + y) n -*(xy) 



Applying Pascal's identity and bringing the leftmost term inside the sum, we obtain 



(a - y ) £(-!)< 



n — i 
i 



(x + y) n - 2i ( xy y 



Finally, we note that for n > 1 the binomial coefficient is when i — n. Therefore, this last expression is 
equal to T(n + 1), and so we are done. □ 

Proof of Theorem^ Begin by expanding equation ijjjj via the binomial theorem to obtain 



Pn(x) = 



ri(-r 2 )^'- fc -r 2 (-r 1 ) 2 "- fc 



On (a;) = 



a 2 "" 


-i 


ri - 




a 2 " 


-l 






Tl - 


'2 


a 2 " 


-1 






T\ - 


''2 


a 2 " 


-1 






n - 


r-2 



k=0 
2 







or 












(T 




V ft, , 





\2 n -k 



2 ,l -X 

E 

fe=0 



(-r 2 ) 2n - fe -(-n) 



2"-fc 



We can now apply lemma to the expressions of the form x n — y n in both P n and Q n . Doing so and 
canceling the factor of r% — r 2 (since b 2 — Aac ^ 0) yields 



Pn(x) 



Qn(x) 



2" -2 

E 



k=0 



2"-fe-2 

E (-i: 

3=0 



;( 2™-fc-j-2 



a 2 "" 1 £ 



fc=0 



2" 



2™-fc-l 

E (-D 

j=0 



2 n - k - j - 1 



(-n-r 2 ) 2n -*- i W- a (r 1 r 4 y +1 
(-ri-^) 2 "-^- 1 ^^ 
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Replacing 7*1 and r^ with their values in terms of the coefficients a, b, c gives the final form 



Pn(x) 



2™-l 2™ 

a x 



Qn(x) = ^ 



k=0 



2™-2 

E 

k=0 

2"-fc-l 



2 n -k-2 

E 

J=0 



2 " fe i 2 \ a k+j b 2"-k-2j-2 ( J+l 



J2 (-1) J ( 2 " ^ J 1 )a k +ib 2n - k -V- 1 ( j 



These equations are the same as equations <|2all and (I2b|) . which was to be proved 



□ 



4 Further Remarks 

Though this paper does not discuss any combinatorial interpretation of the polynomials P n and Q n we 
suspect that their may be some fruitful combinatorial equivalence yet to be discovered. To aid further 
research in this area, we make a couple of observations about P n and Q n . 
Simply iterating Newton's method for the general quadratic gives 



N 



Pn 
Qn 



Pn 



2aP n Q n + bQl Q n+1 



(9) 



where Pq(x) = x and Qq(x) = 1 so that the initial term is zq = x as before. The question is whether the 
numerator and denominator of equation 10 are relatively prime so that we may define P n and Q n recursively 
in the natural way. This is, in fact, true as the following lemma proves. 

Lemma 2. The polynomials P n +\ and Q n +i, defined recursively as 

P n+1 = aPl - cQl 



and 



where Pq{x) 



Q„+i = 2aP n Q n + bQl 
x and Qo(x) — 1, are relatively prime except possibly in the case b 2 — 4ac = 0. 



(For the duration of the statement and proof of Lemma 2, we are suspending the definition of P n and Q n 
given earlier, but it will be an immediate consequence of Lemma 2 that the two definitions agree.) 

Proof. We proceed inductively by assuming that Pi and Qi are relatively prime for all i < n. Assume P n +i 
and Qn+i are not relatively prime to derive a contradiction. Then there exists an irreducible polynomial a 
such that a \ P n +i and a \ Q n +i- If ct \ Q n , then, since a \ aP 2 — cQ^, it follows that a \ P n which contradicts 
the induction hypothesis that P n and Q n are relatively prime. Hence a / Q n . Since a \ 2aP n Q n + bQ^, we 
know that a | 2aP n + bQ n . But then a \ P n + -J^Qn and also a | P n ± y^Q n - Consequently, a divides their 
difference, so a (j^T \/D which only occurs when b 2 — Aac =0. □ 

In combinatorics it is sometimes useful to consider two formal variables x, y which do not commute with 
each other but instead obey yx = qxy where q is another formal variable that commutes with both x and 
y. This approach, due to Schutzenberger 0], is useful in applications such as counting lattice paths. In our 
case, the polynomials P n and Q n can be easily generalized to the non-commuting case. As a generalization 
of the usual binomial coefficient, the g-binomial coefficient is defined as 



and the following g-binomial theorem for non-commuting variables x 1 y due to Schiitzenberger 0] is 

n 

(x + yT = Y, 



k=0 L 



x k y n ~ k 



Analogously to the above, we can then define non-commuting polynomials P' n and Q' n recursively as 

P'n+i = aP n ~ cQn (10a) 



Q'n+l = aP 'nQ'n + aQ'nK + b Qn (10b) 

where Pq — x, Q' = y, and yx = qxy. We conjecture explicit formulas for both P' n and Q n which happen to 
be the same as equations (12 all and (|2b|) except for the presence of a g-binomial coefficient: 

2 n -2 2™-fc-2 

fe=o j=a 

2"— 1 2™-fc-l 

, n(«,y)= E E 
5 Conclusion 

We have taken the initial observation that Newton's method, when applied to quadratics, produces poly- 
nomials with highly composite coefficients and proved an explicit formula for the nth iterate that explains 
this compositeness as a consequence of the inherent compositeness of binomial coefficients. Furthermore, a 
recursive definition and a conjectural non-commutative analogue of the polynomials P n and Q n were noted 
in hopes of spurring further research into finding a combinatorial interpretation. We believe that a proof 
of the non-commutative analogue, as well as the larger issue of finding a combinatorial interpretation, are 
problems which merit further study. It is also worth noting that for higher-degree polynomials, such as 
cubics, no similar phenomena have been found. In particular, the occurrence of coefficients with large prime 
factors indicates that no simple product formulas for the coefficients exist, but this does not rule out the 
existence of more complicated formulas. 
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well as to the computing staff at UW-Madison for providing the computational resources necessary for our 
investigations. 



2 n - k - j - 2 



fc+il,2 n — A— 27— 2 7 + 1 k 2" 

a J o J & x y 



(11a) 



2™ 
k 



-k-j-1 



i fe +^ 2 "- fe -^-VaV 



-k 



(lib) 
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